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DISCUSSION 


Assoc. ASCE.—In this very interesting paper Mr. 
Salvadori has shown that the solution differential equation actually the 
simultaneous solution infinite number algebraic equations and, con- 
versely, for many practical cases, the simultaneous solution finite number 
algebraic equations may accepted lieu the solution differential 
equation. The practical consequences this conclusion are apparent and Mr. 
Salvadori has solved some problems that arise the theory structures. The 
Rayleigh energy method,? widely used for the approximate solution problems 
elastic stability and vibrations, depends the assumption the deflected 
shape the structure. There are some cases (the curved bar, for example) 
which the shape the deflected bar cannot estimated and the solution 
the differential equation required. 

The absence any physical explanation the method finite differences 
somewhat detrimental full appreciation and understanding the method. 
Physically, the method finite differences equivalent substituting less 
complex structure for the ideal mathematical structure. 1920 Heinrich 
proposed that, studying the elastic stability structures, the 
continuous elastic structure replaced one consisting rigid sections 
connected with hinges capable resisting elastic deformation. For example, 
the column Fig. replaced Fig. 13, which parts AB, BC, 


Elastic joint 
Elasticity 


etc., are considered rigid and the elasticity concentrated the hinges. From 
study this model the designer can decide whether the substitute structure 
more rigid less rigid than the original continuous elastic structure, and the 
critical load expected greater than less than that for the continuous 


paper Mario Salvadori was published December, 1949, Proceedings. The number- 
ing footnotes, equations, tables and illustrations this Separate continuation the consecutive 
numbering used the original paper. 

Engr., Bureau Reclamation, Denver, Colo. 
Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
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elastic structure. The designer also interested knowing approximate 
solution gives values the safe side the unsafe side. The problem thus 
longer the realm pure mathematics and enters the domain physics. 

That Mr. Salvadori’s application finite first order differences mathe- 
matically identical with Mr. Hencky’s chain” method apparent from 
comparison the mathematical computations involved. The bar chain 
method for the problem Fig. 11(a) furnishes set linear algebraic equa- 
tions expressing the relation between curvature and moment. Thus, 


and 
2 
The following determinant obtained— 
0..(80) 
which 
The determinantal equations become 


9.37258 


Further have been made torsional equilibrium, curved 
bars, and the elastic equilibrium frames. point emphasized 
that the solution for the substitute structure the 
solution the mathematical continuous structure and only remains for the 
designer decide whether the substitute structure sufficiently similar the 
actual structure for his purposes. From the results given the paper can 
inferred that the bar chain column subjected axial loads less rigid 
structure than the continuously elastic column and that results obtained are 
the safe side. the other hand, for case (Section the bar chain beam 
less rigid more rigid depending whether elastic resistance concentrated 
odd even number hinges. The behavior thus tied with 
the stiffness the substitute structure compared with the ideal mathematical 


structure. 
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Bruno clear presentation the finite differences approach, 
and very logical development valuable method for buckling load calcu- 
lations, are given this paper. One the chief advantages methods such 


this that, although some- 


lengthy, complicated mathe- 
matical manipulation required 
which might overshadow the 


ous steps. fact, once the 
governing differential equation 


9.000 9.373 9.549 9.646 
9.852 9.863 9.867 9.868 


for particular problem established, the buckling load may obtained with 
operations not more complicated than those algebra. 


The method suggested 
useful when possi- 
ble, without too much diffi- 
culty, write differential 
equation (or least 
system differential 
equations) governing the 
problem under consider- 
ation, and the correspond- 
ing boundary conditions. 
some cases, however, 
system equations; ex- 
amples such problems 
are the buckling rein- 
forced plates shells, 
the instability frame- 
works. problems 
some other methods proba- 
bly could more easily 
applied—in particular, 
some the methods de- 
veloped during the investi- 
gations 1947 the 
Polytechnic Institute 
Brooklyn, 
These methods 
are also finite difference 
methods the sense that 


the structure under investigation considered composed several units 


finite dimensions. 


Eng. Specialist, Goodyear Aircraft Corp., Akron, Ohio. 


Methods for the Calculation Elastic Bruno Boley, Journal the 


Aeronautical Sciences, June, 1947, 337. 
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The usefulness the method suggested Mr. Salvadori greatly in- 
creased, course, the extrapolation procedures suggested the paper. 
The examples the paper show that these procedures provide very powerful 
tool. Their chief drawback (as Mr. Salvadori points out) that not 
possible know beforehand whether the value they yield for the buckling 
load higher lower than reality. may mentioned here that the 
methods referred will yield exact answer the analysis 
the units sufficiently rigorous. general, the units may chosen 
ciently simple that accurate analysis possible. assumption 
made (as those methods allow) concerning the deflected shape the structure, 
may shown that the resulting calculated buckling load will too high. 

obtain better idea the accuracy afforded the extrapolation pro- 
cedures, necessary that additional problems solved involving different 
types structures. any one problem, however, good estimate the 
correct buckling load may probably obtained some graphical method 
extrapolation, such the rather obvious one (for the case the Euler buckling 
straight column) outlined herein. The values Table were given 


divided into equal segments, and the value obtained the 
extrapolation—that is, 


for the case which the bar assumed 


becomes very large, the values and approach the correct value, 


Thus, plots are made against against and against 
the common intersection the three resulting curves will represent the 
correct value, Furthermore, this intersection will occur line which 
the abscissa equal the ordinate. These plots are given Fig. for the 
present example. 

number such schemes may devised either improve the accuracy 
the method estimate the error introduced considering the structure 
comprising finite number points. many cases, however, the extrapola- 
tion procedure leads such good results that improvement may required 
for engineering accuracy. 


Jun. ASCE.—The use difference equations ob- 
tain approximate solutions boundary value problems science and engineer- 
ing not new;* but the author’s clear exposition certain make the method 
familiar engineers. 

Prof. Civ. Eng., Annamalai Univ., Annamalainagar, South India. 
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The writer has undertaken check the validity the extrapolation pro- 
cedures discussed the author and study the effect varying boundary 
conditions convergence. The example chosen the classical Euler problem, 
the strut being hinged one end and re- 


strained position and direction the TABLE 
other end. The exact for CoMPUTATIONS 
this case 20.16 Table 
order differences. 


with Example Mr. Salvadori’s 

paper reveals that convergence very 

much slower for the boundary conditions assumed the writer. However, 
h?-extrapolations give surprisingly accurate results even with very few trials. 


power the Richardson method” extrapolation 
the computation buckling loads demonstrated this paper. Although 
the solution elastic stability problems finite differences has been subject 
much study many eminent workers, without the use the method 
extrapolation, the computations become very tedious when accurate solution 
required. applying the method extrapolation, Mr. Richardson as- 
sumed that the error the buckling load also the h?-type which the 
author has argued not necessarily true. However, solving many nu- 
merical examples, Mr. Salvadori has shown that the assumption did work; 


but, because the lack proof, had limit his conclusion adding the 


phrase experimental Nevertheless, the fact that the 
assumption did work out well all these cases appears more than mere 
coincidence. 

Mr. Richardson’s assumption can actually proved rigorous 
least for the cases where the moment inertia constant. Consider the case 
the buckling column with hinged ends under axial compression. The 
differential equation 


responding difference problem 


and the approximate value obtained solving the difference equation 
(Eq. 85) rather than the differential equation (Eq. 84). 


Elastic Timoshenko, McGraw-Hill Book Co., Inc., New York, Y., 
1st Ed., 1936, 89. 

Associate Prof. Aeronautical Eng., Daniel Guggenheim School Aeronautics, College Eng., 
New York Univ., New York, 
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The finite difference equation (Eq. 85) can solved assuming the 
solution the form: 


Substituting Eq. Eq. 85, 
The form Eq. requires that 
Three forms the solution for may obtained depending whether 


can shown that the resulting solution for cannot satisfy the boundary 
conditions and therefore that these cases will not considered. When 


therefore, 


Substituting values from Eq. Eq. 86, the values can then 
determined from the boundary conditions. this procedure proves rather 
laborious, the solution simplified first determining the value from 
the boundary conditions and then solving from Since can either 
positive negative, the general solution for 


satisfy the boundary conditions, necessary that 


which any integer, and the instant buckling. Since 
tained from Eq. 94; thus: 


— 


— 
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2 


2 
results from the fact that thus proving that 


the may noted that for columns with other boundary condi- 
tions, the differential equation the form given Eq. 84, Eq. 
shows that always equal and therefore the proof holds equally true 
these cases. general, the governing differential equations are those 
with constant coefficients, the corresponding finite difference equations can 
solved the same symbolic method the case the differential equations 
and the final solutions for have the same form those for 
will always equal and the relation given Eq. can always 
proved correct. 


consideration stability problems different 
fields engineering steadily becoming greater importance. The value 
this paper lies two facts: (1) The finite differences method (which had 
not been thoroughly investigated) presented very broad and clear 
manner; and (2) through his simple extrapolation formulas, the author increases 
the power the method large extent. 

The worker engaged the investigation characteristic value problems 
will thankful the author for considering many aspects the method 
(such different coordinate systems, different order operators, procedures 
for numerical computations, types asymptotic approach, etc.), and will 
notice that the method does not require any guessing the deflection line (or 
surface), such required the energy method and other approximate 
Since complicated boundary conditions (for example, elastically 
supported can expressed finite differences, the method ap- 
plicable larger variety problems than indicated the paper. The 
extrapolation formulas show amazing power, can seen the many 
illustrative examples. their use the numerical computations are made 
less cumbersome, with the saving time and energy. 


Instr. Civ. Eng., Columbia Univ., New York, 
Internationale Vereinigung fiir und Hochbau, Bd. 1947, 237. 
der rechteckigen Platte bei veriinderlicher Shizuo Ban, 
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There are, however, many problems elastic stability which the energy 
method leads good results shorter time than any other approach. For 
the plate Fig. exposed thrust triangular shape 
along the simply supported edges The load-free edges, are built and 
free, respectively. Apply the energy method, assuming sinusoidal deflection 
function the direction the thrust and deflection function the type— 


the direction perpendicular the thrust. the ratio a/b 1.635, Eq. 
yields 1.636 the first approximation and practically the same value 
the second approximation, the difference being noticeable only the fifth 
digit. Using the finite differences method with first order operators, 1.14 
1.42; the error still and find was necessary solve 
fourth degree equation. 

Although this example not typical, far convergence concerned, 
general, the investigator dealing with stability 
plates, interested finding the ratio a/b the 
edges the plate for which minimum. The 
buckling equation being established, is, general, 
easier find (minimum) the energy method, 
because the ratio a/b appears only the diagonal 
members the determinant, whereas, using finite 
differences, this ratio appears most elements the 
determinant. 

The extrapolation formulas can applied 
solving not only the characteristic value problems, 
but also the boundary (or initial) value problems. 
For given problem the type approach depends 
the choice the type coordinates. Thus the 

problem skew plate with equal sides under 

uniform load and simply supported along all four 

edges, the true value for the deflection can squeezed, using different types 
coordinates. 

Using skew coordinates for the analysis 60° skew plate, 
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Beulung von auf einseitigen, Druck beanspruchten Platten,” 
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for being the flexural rigidity the Applying the extra- 
polation formula probably better value 0.0396 results. 

Using triangular coordinates, computed the values 


rounded off here three figures. this case the extrapolation formula (Eq. 


4 
41) gives 0.0408 The average the two extrapolated values 


4 
0.0402 and, plausibly, does not differ more than approximately from 


the true value the deflection the center. This example illustrates, the 
same time, how important choose appropriate type coordinates, 
emphasized the author. 

Extremely good approximations are obtained, using extrapolation formulas, 
computing the deflections uniformly loaded square plates. the case 


4 4 
simple supports," 0.06250 and 0.06445 The extrapolated 


4 
value (using Eq. 41) equals 0.06510 and differs only 0.2% from 


4 
the exact value 0.06496 the case built-in 


4 4 
0.00180 and 0.00143 spite the fact that these values 


4 
differ from the exact value 0.00128 41% and 12%, respectively, 


4 
the extrapolation formula (Eq. 41) leads 0.00131 which only 


2.3% error. 

the case transverse concentrated loads applied plate, the extra- 
polation formulas given the author cannot used, because the error 
due not only the substitution finite differences for differentials, but also 
the spreading the load over certain area.* investigating skew plates 
under central concentrated the writer had use higher order differ- 
ences obtain satisfactory approximations. might worthwhile investi- 
gate the possibility developing extrapolation formulas problems this type. 


Mario Assoc. ASCE.—Mr. Silverman points out 
that the Hencky “bar chain” method equivalent the finite differences 
approach and emphasizes the advantages physical picture solving 
engineering problems. The choice between physical mathematical 


Theorie elastischer Gewebe und ihre Anwendung auf die Berechnung biegsamer Platten,’’ 
Marcus, Springer, Berlin, Vol. Ed., 1932, 61. 

Untersuchungen die Durchbiegung Platten unter 
Polytechnique Féderale, Verlag Leemann, Ziirich, 1950. 
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approach rests essentially the mental attitude the worker and the 
relative efficiency the two types approach. this connection will 
suffice note that the physical approach Mr. Hencky would never lead 
the investigation the error and the discovery the extrapolations, 
without which his method excessively laborious and inefficient. The prob- 
lems dealt with the paper illustrate the decided superiority the mathe- 
matical approach the present case. other cases, combination physical 
intuition and mathematical knowledge usually gives the best results. 

The graphical “extrapolation the extrapolations,” suggested Mr. 
Boley, although without sound mathematical justification, may prove useful. 
can applied any converging series approximate values and has the 
advantage bringing the asymptotic value the approximations point 
45° line. doubtful whether the graphical extrapolation would give 
accurate results the case oscillating approximations. 

Mr. Ramaswamy’s example, which checks unpublished result obtained 
the writer, clearly indicates the powerfulness the extrapolation procedure. 
interesting add that the (h? applied the approxi- 
mations gives (in this case) 20.166, with error 0.03%. 

The writer grateful Mr. Herrmann for showing the applicability 
the extrapolations problems other than those treated the paper. The 
extrapolations can actually used connection with any method successive 
approximations which: (a) The successive values the unknown depend 
increasing integer and (b) the error expandable into series negative 
even powers Various procedures this type, which the extrapolations 
can usefully applied have been discussed elsewhere the 

The proof, given Mr. Wang, that the error the characteristic value 
Euler’s problem the h?-type represents important contribution 
the theory the finite differences method. opens the path more 
general proof which valid for all ordinary equations with constant coefficients. 

The writer unable follow the reasoning the last paragraph Mr. 
Wang’s discussion, however, and does not see that proves Eq. 
correct for all equations with constant coefficients. fact, Eq. was derived 
for the particular case Eq. and this equation does not hold even for all 
beams with constant moment inertia. The following fundamental theorem 
may instead proved, the error-formula (Eq. 34) given the writer: 


Theorem. The error the characteristic value linear ordinary 
differential equations with constant coefficients evaluated finite differ- 
ences the h?-type. 


For this purpose notice that, according the theorem Richardson, 

and and hence are the Therefore, the condition: 
sufficient condition (although not necessary condition) for the 
h?-type; but the solution ordinary difference problem with constant 


January, 1949. 
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coefficients identical with solution the corresponding differential prob- 
and, hence, Eq. (the condition Eq. 100), satisfied when 
and have constant coefficients. Thus, for equations with constant coeffi- 
cients the error reduces the form: 


now obvious that the extrapolations must give good results con- 
nection with constant coefficient equations and that the accuracy the 
extrapolations must improve taking into account larger number terms 
the series for the error. For example, the writer has found that, large 
number cases, the (h? h‘)-extrapolation gives better results than the 
h?-extrapolation connection with constant coefficient equations. 

also interesting note that the behavior the successive approxi- 
mations monotonic all the constant coefficient examples solved the 
writer, whereas the behavior seldom monotonic for equations with variable 
coefficients. For this last type equation the sufficient condition Eq. 100, 
general, not satisfied; and this fact may proved checking the values 
the approximate solution the pivotal points, obtained finite differences, 
note that they differ from the values the integral the differential equa- 
tion the same points. For example, the differential problem for the buckling 
strut, built one end and free the other, under the action its 
own weight may written the form: 


and 
and has series solution: 


(S. Timoshenko? develops Eq. 102; here dy/dz,z q/(E1).) 
dividing the strut into two equal parts, the ratio the middle point, 
slope (L/2), the end slope, (0), becomes (with 7.81 and, 
hence, 0.977): 0.842/1 0.842. The corresponding dif- 
erence problem given the equations: 


and 


(107) 


und partiellen Differenzengleichungen der Baustatik,” Fr. Bleich and 
Melan, Springer, Berlin, 1927. 


For Eq. 105 and Eq. 106 give rise the system two equations: 
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and 


273 


and, with this value for the ratio becomes 


Hence, the approximate differs from the true and Eq. 100 not necessarily 
satisfied. 

Finally, must noted that, for equations with variable coefficients, 
increase the number terms the error series taken into account usually 
reduces the accuracy the extrapolations obtained means all the succes- 
sive approximations. Therefore, these results seem indicate that the error 
the characteristic values equations with variable coefficients not the 
h?-type. 

The reason why the extrapolations work well connection with any 
monotonic series approximations remains investigated, although this 
result could well purely empirical character. 
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